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Abstract 

BACKGROUND 

There is no human epidemiology data available to estimate the heavy ion 
cancer risks experienced by astronauts in space. Studies of tumor induction in 
mice are a necessary step to estimate risks to astronauts. Previous 
experimental data can be better utilized to model dose response for heavy ion 
tumorigenesis and plan future low dose studies. 

DOSE RESPONSE MODELS 

The Harderian Gland data of Alpen et al.[l-3] was re-analyzed [4] using non- 
linear least square regression. The data set measured the induction of 
Harderian gland tumors in mice by high-energy protons, helium, neon, iron, 
niobium and lanthanum with LET's ranging from 0.4 to 950 keV/micron. We 
were able to strengthen the individual ion models by combining data for all ions 
into a model that relates both radiation dose and LET for the ion to tumor 
prevalence. We compared models based on Targeted Effects (TE) to one 
motivated by Non-targeted Effects (NTE) that included a bystander term that 
increased tumor induction at low doses non-linearly. When comparing fitted 
models to the experimental data, we considered the adjusted R 2 , the Akaike 
Information Criteria (AIC), and the Bayesian Information Criteria (BIC) to test for 
Goodness of fit. In the adjusted R 2 test, the model with the highest R 2 values 
provides a better fit to the available data. In the AIC and BIC tests, the model 
with the smaller values of the summary value provides the better fit. The non- 
linear NTE models fit the combined data better than the TE models that are 
linear at low doses. We evaluated the differences in the relative biological 
effectiveness (RBE) and found the NTE model provides a higher RBE at low dose 
compared to the TE model. 

POWER ANALYSIS 

The final NTE model estimates were used to simulate example data to consider 
the design of new experiments to detect NTE at low dose for validation. Power 
and sample sizes were calculated for a variety of radiation qualities including 
some not considered in the Harderian Gland data set and with different 
background tumor incidences. We considered different experimental designs 
with varying number of doses and varying low doses dependant on the LET of 
the radiation. The optimal design to detect a NTE for an individual ion had 4 
doses equally spaced below a maximal dose where "bending" due to cell 
sterilization was < 2%. For example at 100 keV/micron we would irradiate at 
0.03 Gy, 0.065 Gy, 0.13 Gy, and 0.26 Gy and require 850 mice including a control 
dose for a sensitivity to detect NTE with 80% power. Sample sizes could be 
improved by combining ions similar to the methods used with the Harderian 
Gland data. 

Color Categories 

Uncertainty Reduction/Risk Mitigation Category Color = Green 
Mechanistic/Descriptive Category Color = 

Blue Bullets Progress 

• Applying a NTE model motivated by the dose response observed for 
bystander effects and genomic instability in cell culture, we show that the 
NTE model provides a superior fit to the dose response for tumors in mice 
based on several model ranking tests. These results add important empirical 
evidence in support of the NTE model based on in vivo data for tumor 
responses. 

• We were able to fit models for the LET dependence of tumors over a broad 
range from protons to heavy ions and describe a dose dependent RBE 
applicable to low dose exposures. 

• Power analysis based on the Harderian gland data suggest the optimal 
number of mice to be studied for future dose response experiments and 
suggest sample size reductions will occur when results for several radiation 
qualities are combined into a single model. 


Red Bullets Gaps in Progress and Knowledge 

• The relative contribution to cancer risks from targeted effects and non- 
targeted effects remains elusive with too few experiments designed to test 
the shape of the dose response at low doses (<0.3 Gy) applicable to space 
missions. 

• Only a few murine model tumors have been studied with only a few ion 
types. The paucity of data limits the building and testing of models of cancer 
risk from space radiation. 

• The usage of human cell culture models in 2D or 3D is needed to support the 
applicability of murine models to human risk prediction, however much work 
remains in making the necessary connections. Of importance is the need for 
more expansive data set on radiation quality at a variety of low doses to 
understand the shape of the dose response for cancer processes induced by 
heavy ions. 


Description of Alpen data. 


Ion 

Energy 

(MeV/u) 

LET 

(keV/nm) 

Hydrogen 

250 

0.4 

Helium 

228 

1.6 

Neon 

670 

25 

Iron 

600 

193 

Iron 

350 

253 

Niobium 

600 

464 

Lanthanum 

593 

953 


Measured Harderian gland tumors in female BGCF^Anl mice. 

A variety of heavy ions were accelerated at the Lawrence Berkeley 
National Lab for the exposures. 

Mice were irradiated at 100 to 120 days. 

Tumor prevalence was determined 16 months after irradiation. 

Tumor appearance was accelerated using pituitary implants 
o Reduced the effects of competing tumor risks and reduced costs of 
experiments (sacrifice at 600 d compared to lifespan of ~900 d) 

Fe nuclei showed same tumor response with or without the 
pituitary implants. 

y-rays were used as a low-LET reference radiation 


y-ray TE model relating dose to tumor prevalence 


Applications of Models for Risk Estimation 



Cross-over dose is defined as 
the dose where TE = NTE 



A dose dependent relative biological 
effectiveness (RBE) is found in the 

NTE model 


P u = P, +[aD + pD 2 ]e XD 


Dose induction Cell sterilization 

term term 



Estimate 

Standard 

Error 

Po 

2.64 

0.77 

a 

4.31 

3.09 

3 

5.08 

2.25 

A 

0.264 

0.055 



Heavy ion TE model relating dose and LET to tumor prevalence 


P TE = P 0 +[a(L)D]e~ ML)D 

a(L) -a 0 + a x L exp(-a 2 Z) 

A (L) - A 0 + AyL Qxp(-A 2 L) 



Estimate 

Standard 

Error 

Po 

2.8 

0.48 

a o 

12.39 

2.2 

a i 

1.17 

0.21 

a 2 

0.00369 

0.00064 

K 

0.0504 

0.0298 

K 

0.0068 

0.0029 


0.0033 

0.002 
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Adjusted r 2 = 0.9296, AIC = 210, BIC = 222.1 

Heavy ion NTE model relating dose and LET to tumor prevalence 

P N ^=P 0 +[a(L)D]e- iWD 

+ [K t L exp(-x: 2 Z,)] 0 (Z), s )e~ ML)D 
©( D lh ) represents a step function 
with threshold dose ( D th = 0.01 Gy) 
for the bystander effect. 
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Estimate 

Standard Error 

Po 

2.44 

0.41 

a o 

13.60 

1.87 

a l 

0.64 

0.21 

a 2 

0.00353 

0.00081 
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0.055 
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0.0055 
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Adjusted r 2 = 0.9514, AIC = 195.9, BIC = 211.6 


P-values are 0.056 and 0.03 respectively. 


Power Analysis Specifications 


• The study focuses on low doses where the cell sterilization 
term can be ignored. 

• Parameter estimates from the NTE model were used to 
generate simulations for each ion individually. 

• Simulations were analyzed using generalized linear models 
with binomial errors following the low dose model. 

• Power was determined by the ability to detect if NT Effects 
(/c) are significant at low dose 

• Doses were chosen at OGy and 3 or 4 irradiation doses 
chosen as seen in the figure with sample sizes such that all 
doses have equal binomial variances. 

• Dosing Schemes for irradiated mice: 

1. Three radiation doses (Dl, D2, and D3 from figure) were 
considered and d ~ 2%. 

2. Three radiation doses were considered and d ~ 1%. 

3. Four radiation doses (Dl, D2, D3, and D4 from figure) 
were considered and d ~ 2%. 


Equation: The individual NTE Model 

P,„=P o +[aD + K 0(DJ]e- AD 

Equation: The linear low dose Model 

P ion =P Q +ccD + K®(D th ) 

Table: The LET values chosen for the 
study and parameter estimates 


from the NTE model 


a 

A 

K 

LET = 70 

48.59 

0.317 

5.47 

LET = 100 

58.56 

0.372 

6.82 

Fe, 

76.01 

0.449 

8.62 

LET = 193 





Figure: Representation of how 
doses were chosen. Each additional 
dose is half way between the 
previous dose and 0. 


Table: Sensitivity of sample size to detect NTE with 80% power to the scheme 
choice and background %-Prevalence (P 0 ). 



Scheme 1 

Scheme 2 

o/- 

Scheme 3 

Scheme 3 

Scheme 3 


P 0 = 2.44 

P 0 = 2.44 

P 0 = 2.44 

Po = 5 

o 

rl 

II 

o 

Q_ 

Fe, 

900 

1031 

690 

810 

1069 

LET = 193 






LET = 100 

1200 

1180 

850 

1061 

1810 

LET = 70 

1620 

1620 

1101 

1600 

>2100* 

Total** 

3120 

3231 

2041 

2871 

>4379 


*Power = 68.27% **Assuming controls same across the 3 ions 



Dose 

"d" is %-difference between linear 
low dose model (light green) and the 
individual NTE models (dark green) 


Power Analysis Conclusions 

• Dosing Scheme 3 was optimal. 

• More mice are needed as the value of P 0 increases. 

• The sample size estimates to test for NTE for single 
ions will over-estimate sample size requirements 
when results from several ions are combined into a 
model that describes the LET dependence of the 
tumor response. 


Figure: Doses and sample sizes for Scheme 3 
by individual ions 
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Source of tumor data for heavy ions 

[1] Fry RJM, Powers-Risius P, Alpen EL, Ainsworth EJ. High LET radiation carcinogenesis. Radiat. Res. 104: S188-195 (1985). 

[2] Alpen EL, Powers-Risius P, Curtis SB, DeGuzman R: Tumorigenic potential of high-Z, high-LET charged particle radiations. Radiot. Res. 88: 132-143 

(1993). 

[3] Alpen EL, Powers-Risius P, Curtis SB, DeGuzman R, Fry RJM. Fluence-based relative biological effectiveness for charged particle carcinogenesis in 

mouse Harderian gland. Adv. Space Res. 14: 573-581 (1994). 

Related Literature 

[4] Cucinotta FA, and Chappell LJ, Non-targeted Effects and the Dose Response for Heavy Ion Tumor Induction. Mut. Res. 687: 49-53 (2010). 

[5] Weil MM, Bedford JS, Bielefeldt-Ohmann H, Ray AF, Gernik PC, Ehrhart EJ, Falgren CM, Hailu F, Battaglia CLR, Charles C, Callan MA, and Ullrich RL. 

Incidence of acute myeloid leukemia and hepatocellular carcinoma in mice irradiated with 1 GeV/nucleon 56 Fe ions. Radiat. Res. 172: 213-219 
(2009). 

[6] Dicello JF, Cucinotta, FA, Gridley, D, Mann J, Moyers MF, Novak GR, Piantadosi S, Ricart-Arbona R, Simonson DM, Williams, JR, Zhang Y, Zhou H, and 

Huso D In vivo mammary tumorigenesis in the sprague-dawley rat and microdosimetric correlates. Phys. Med. Biol. 49: 3817-3830 (2004). 

[7] Cucinotta FA and Wilson JW. Initiation-promotion model for tumor prevalence from high energy and charge radiation. Phys. Med. Biol. 39: 1811- 

1831 (1994). 



















